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Abstract. It has been known for a long time that the cross section for bottom-quark b production 
at hadron collider energies exceeds theoretical expectations. An additional contribution from pair- 
production of light gluinos g, of mass 12 to 16 GeV, with two-body decays into bottom quarks 
and light bottom squarks b, helps to obtain a b production rate in better agreement with data. The 
masses of the g and b are restricted further by the ratio of like-sign to opposite-sign leptons at hadron 
colliders. Constraints on this scenario from other data are examined, and predictions are made for 
various processes such as Upsilon decay into Fs. 


INTRODUCTION 

The cross section for bottom-quark production at hadron collider energies exceeds the 
central value of predictions of next-to-leading order (NLO) perturbative quantum chro¬ 
modynamics (QCD) by about a factor of two [|I], ^ . This longstanding discrepancy has 
resisted fully satisfactory resolution within the standard model [|^]. The NLO contribu¬ 
tions are large, and it is not excluded that a combination of further higher-order effects 
in production and/or fragmentation may resolve the discrepancy. However, the disagree¬ 
ment is surprising because the relatively large mass of the bottom quark sets a hard 
scattering scale at which fixed-order perturbative QCD computations of other processes 
are generally successful. The photoproduction cross section at DESY’s HERA [Q] and 
and the cross section in photon-photon reactions at CERN’s LEP [Q also exceed NLO 
expectations. The data invite the possibility of a contribution from “new physics". 


SUPERSYMMETRY INTERPRETATION 

The properties of new particles that can contribute significantly to the bottom quark 
cross section are fairly well circumscribed. To be produced with enough cross section the 
particles must interact strongly and have relatively low mass. They must either decay into 
b quarks or be close imitators of h’s in a variety of channels of observation. They must 
evade constraints based on precise data from measurements of decays at CERN’s LEP 
and SLAC’s SEC, and from many lower-energy collider experiments. The minimal 

supersymmetric standard model (MSSM) is a favorite candidate in many quarters for 
physics beyond the standard model. It offers a well-motivated theoretical framework 
and is reasonably well-explored phenomenologically. An explanation within the context 
of the MSSM [0 can satisfy all of the stated criteria. 





In Ref. [|p, the existenee is assumed of a relatively light eolor-oetet gluino g (mass 
~ 12 to 16 GeV) that deeays with 100% branehing fraetion into a bottom quark b and 
a light eolor-triplet bottom squark b (mass ~ 2 to 5.5 GeV). The g and the b are the 
spin-1/2 and spin-0 supersymmetrie partners of the gluon (g) and bottom quark. In this 
seenario the b is the lightest SUSY partiele, and the masses of all other SUSY partieles 
are arbitrarily heavy, i.e., of order the eleetroweak seale or greater. The b is either long- 
lived or decays via R-parity violation into a pair of hadronic jets. Improved agreement is 
obtained with hadron collider rates of bottom-quark production, and several predictions 
are made that can be tested readily with forthcoming data. 


Differential Cross Section 

The light gluinos are produced in pairs via standard QCD subprocesses, dominantly 
g-|-g—i>g-|-g at Tevatron and Large Hadron Collider (LHC) energies. The g has a 
strong color coupling to Vs and Vs and, as long as its mass satisfies nig > nih + mj^, the 
g decays promptly to b + b. The magnitude of the b cross section, the shape of the Vs 
transverse momentum prb distribution, and the CDF measurement of 5^ — 5® mixing 
are three features of the data that help to establish the preferred masses of the g and b. 
In Ref. [||], contributions are included from both q + q^ g + g and g + g —g + g. The 
subprocess g + b ^ g + b contributes insignificantly. Shown in Fig. 1 is the integrated 
Prt> distribution of the b quarks that results from g ^ b + b, for nig =14 GeV and = 
3.5 GeV. The results are compared with the cross section obtained from next-to-leading 
order (NLO) perturbative QCD and CTEQ4M parton distribution functions (PDF’s) [|^, 

with mt, = 4.75 GeV, and a renormalization and factorization scale p = \Jfn^+p\i,- 

SUSY-QCD corrections to bb production are not included as they are not available and 
are generally expected to be somewhat smaller than the standard QCD corrections. A 
fully differential NLO calculation of g-pair production and decay does not exist either. 
Therefore, the g-pair cross section is computed from the leading order (LO) matrix 

element with NLO PDF’s [0, = ^nt^ + pj^, and a two-loop expression for the strong 

coupling a^. To account for NLO effects, this g-pair cross section is multiplied by 1.9, 
the ratio of inclusive NLO to LO cross sections 0. 

A relatively light gluino is necessary in order to obtain a bottom-quark cross section 
comparable in magnitude to the pure QCD component. Values of mg ~ 12 to 16 GeV are 
chosen because the resulting g decays produce prb spectra that are enhanced primarily 
in the neighborhood of ~ nig where the data show the most prominent enhancement 
above the QCD expectation. Larger values of nig yield too little cross section to be of 
interest, and smaller values produce more cross section than seems tolerated by the ratio 
of like-sign to opposite-sign leptons from b decay, as discussed below. The choice of 
has an impact on the kinematics of the b. After selections on large values of 
reduce the cross section and, in addition, lead to shapes of the prb distribution that agree 
less well with the data. 

After the contributions of the NLO QCD and SUSY components are added, the mag¬ 
nitude of the bottom-quark cross section and the shape of the integrated p™™ distribution 
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FIGURE 1. Bottom-quark cross section in pp collisions at y/S =1.8 TeV for pjb > with a gluino 

of mass ifig = 14 GeV and a bottom squark of mass = 3.5 GeV. The dashed curve is the central value 
of the NLO QCD prediction. The dotted curve shows the pr spectrum of the b from the SUSY processes. 
The solid curve is the sum of the QCD and SUSY components. The shaded band represents an uncertainty 
of roughly ±30% associated with variations of the renormalization and factorization scales, the b mass, 
and the parton densities. The rapidity cut on the h’s is |y*| < 1. Data are from Ref. 1. 


are described well. Very good agreement is obtained also with data from the UAl exper¬ 
iment (not shown) [@]. The SUSY process produces bottom quarks in a four-body final 
state and thus their momentum correlations are different from those of QCD. Angular 
correlations between muons that arise from decays of b"s have been measured [0, 

The angular correlations between h’s in the SUSY case are nearly indistinguishable from 
those of QCD once experimental cuts are applied. 

The energy dependence of the bottom cross section is a potentially important con¬ 
straint on models in which new physics is invoked to interpret the observed excess bot¬ 
tom quark yield. Since the assumed g mass is larger than the mass of the b, the g pair 
process will turn on more slowly with energy than pure QCD production of bb pairs. 
The new physics contribution will depress the ratio of cross sections at 630 GeV and 
1.8 TeV from the pure QCD expectation. An explicit calculation with CTEQ4M parton 
densities and the b rapidity selection |y| < 1, yields a pure QCD prediction at NLO of 
0.17 +/- 0.02 for = 10.5 GeV, and 0.16 +/- 0.02 after inclusion of the gluino pair 
contribution. Either of these numbers is consistent with forthcoming data from CDE on 






this ratio [jH)]. 


Same-sign to Opposite-sign Leptons 


Since the f is a Majorana particle, its decay yields both quarks and antiquarks. 
Gluino pair produetion and subsequent deeay to Vs will generate bb and bb pairs, as 
well as the bb final states that appear in QCD produetion. When a gluino is highly 
relativistie, its helieity is nearly the same as its ehirality. Therefore, seleetion of g’s 
whose transverse momentum is greater than their mass will reduce the number of like- 
sign Vs. In the intermediate pr region, however, the like-sign suppression is redueed. 
The euts ehosen in eurrent hadron eollider experiments for measurement of the ratio of 
like-sign to opposite-sign muons result in primarily unpolarized g’s, and, independent 
of the b mixing angle, an equal number of like-sign and opposite-sign Vs is expeeted at 
production. The SUSY mechanism leads therefore to an inerease of like-sign leptons in 
the final state after semi-leptonic deeays of the b and b quarks. This inerease eould be 
eonfused with an enhaneed rate of mixing. 

Time-integrated mixing analyses of lepton pairs observed at hadron colliders are 
interpreted in terms of the quantity % = fdXd + fsXs, where fd and fs are the fraetions 
of Bj and 5® hadrons in the sample of semi-leptonie B deeays, and %/ is the time- 
integrated mixing probability for 5®. Conventional bb pair produetion determines the 

quantity LSc — 2%(1 — %), the fraction of bb pairs that decay into like-sign leptons. The 
SUSY mechanism leads to a new expression 
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The faetor 1/2 arises because N{bb -\-bb) ~ N{bb) in the SUSY mechanism for the 
seleetions on pxb made in the CDF run I analysis. Defining G = Ogg/Oqcd, the ratio of 
SUSY and QCD bottom-quark eross seetions after euts, and solving for the effeetive 
mixing parameter, one obtains 
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The CDF measurement is interpreted in Ref. [Q] as a determination of %eff = 0.131 ± 
0.02±0.016 [0]. This number is marginally larger than the world average value % = 
0.118 ± 0.005 [|T^], assumed to be the contribution from the pure QCD component only. 

The ratio G is determined in the region of phase spaee where the measurement is 
made 0, with both final Vs having pjb > 6.5 GeV and rapidity \yb\ < 1. With = 
3.5 GeV, G = 0.37 and 0.28 for gluino masses nig = 14 and 16 GeV, respeetively. The 
predictions are Xeff = 0.17 ± 0.02 for nig = 14 GeV, and Xeff = 0.16 ± 0.02 with nig = 
16 GeV. Additional theoretical uncertainties arise beeause there is no fully differential 
NLO ealeulation of gluino produetion and subsequent deeay to Vs. The ehoiee mg > 12 
GeV leads to a ealeulated Xeff that is eonsistent with the data within experimental and 
theoretieal uneertainties. With Ogg/Oqcd ~ 1/3, the mixing data and the magnitude and 
PT dependenee of the b production cross section can be satisfied. 










CONSTRAINTS FROM OTHER DATA 


An early study by the UAl Collaboration [ p^ exeludes g's in the mass range 4 < nig < 
53 GeV, but it starts from the assumption that there is a light neutralino whose mass is 
less than the mass of the gluino. The eonelusion is based on the absenee of the expeeted 
deeay g^q + q+ where represents the missing energy assoeiated with the x®. In 

the seenario diseussed above, this deeay proeess does not oeeur sinee the bottom squark 
is the LSP, the SUSY partiele with lowest mass, and the mass is presumed to be large 
(z.e., > 50 GeV). An analysis of 2- and 4-jet events by the ALEPH eollaboration [ |T4| ] 
disfavors g’s with mass nig < 6.3 GeV but not g’s in the mass range relevant for the 
SUSY interpretation of the bottom quark produetion eross seetion. A similar analysis 
is reported by the OPAL eollaboration [|T3|]. A light b is not exeluded by the ALEPH 
analysis. The exelusion by the CLEO eollaboration [ |T^ ] of a h with mass 3.5 to 4.5 GeV 
does not apply sinee their analysis foeuses only on the deeays b civ and b —> ci. 
The b need not deeay leptonieally nor into eharm. On the other hand, these data might 
be reinterpreted in terms of a bound on the R-parity violating lepton-number violating 
deeay of b into c£. It would be interesting to study the hadronie deeays b —> cq, with 
q = d or s, and b ^ us with the CLEO data. The DELPHI eollaboration ’s0 seareh 
for long-lived squarks in their yy event sample is not sensitive to < 15 GeV. 

There are important eonstraints on eouplings of the bottom squarks from preeise 
measurements of ZP deeays. A light b would be ruled out unless its eoupling to the tP is 
very small. The squark eouplings to the ZP depend on the mixing angle 0/,. As deseribed 
in Ref. [jT^], if the light bottom squark (hi) is an appropriate mixture of left-handed 
and right-handed bottom squarks, its lowest-order (tree-level) eoupling to the ZP ean be 
arranged to be small if sin^Gz, ~ 1/6. The eouplings and survive, where ^2 
is the heavier bottom squark. However, as long as the eombination of the masses -|- 


m 


bi 


is less than the maximum eenter-of-mass energy explored at LLP, these eouplings 
present no diffieulty. This eondition roughly implies > 200 GeV. However, mueh 

lower masses of ^2 might be tolerated. A eareful phenomenologieal analysis is needed of 
expeeted ^2 deeay signatures, along with an understanding of deteetion effieieneies and 
expeeted event rates, before one knows the admissible range of masses eonsistent with 
LLP data. At higher-order, unless the ^2 mass is of order 100 GeV, eontributions from 
loop proeesses in whieh light gluinos are exehanged may produee signifieant deviations 
from measurements of the ratios At, the forward-baekward b asymmetry at the Z®, and 
Rh, the hadronie branehing ratio of the ZP into b quarks [[1^, |20| ]. 

Bottom squarks make a small eontribution to the inelusive eross seetion for e+e —> 
hadrons, in eomparison to the eontributions from quark produetion, and bb resonanees 
are likely to be impossible to extraet from baekgrounds [ ^Tj ] . The angular distribution of 
hadronie jets produeed in e+e annihilation ean be examined in order to bound the 
eontribution of sealar-quark produetion. Spin-1/2 quarks and spin-0 squarks emerge 
with different distributions, (1 ieos^O), respeetively. The angular distribution measured 
by the CELLO eollaboration [|2^1 is eonsistent with the produetion of a single pair of 
oharge-1/3 squarks along with five flavors of quark-antiquark pairs. A new examination 
of the angular distribution with greater statisties would be valuable. 




PREDICTIONS AND IMPLICATIONS 


T Decay into Bottom Squarks 


If the bottom squark mass is less than half the mass of one of the Upsilon states, 
then Upsilon decay to a pair of bottom squarks might proceed with sufficient rate 
for experimental observation or exclusion of a light bottom squark. In Ref. p3|], the 
expected rate for T —> bb* is computed as a function of the masses of the bottom 
squark and the gluino. The mass of the gluino enters because the g is exchanged in 
the decay subprocesses. The electronic width of the T is used as the source of absolute 
normalization. 

The data sample is largest at the T(4S). For a fixed gluino mass of 14 GeV, the 
branching fraction into a pair of bottom squarks is about 10 for = 2.5 GeV, and 
about 10^^ for = 4.85 GeV. A sample as large as 10,000 may be available in current 
data from runs of the CLEG detector. 

The predicted decay rates for the T{nS), n = 1, 3 are shown in Fig. 2 These curves 
can be read as predictions of the width for the corresponding values of and ntg, or 
as lower limits on the sparticle masses given known bounds on the branching fractions. 
The current experimental uncertainties on the hadronic widths of the T’s are compatible 
with the range of values of and nig favored in the work on the bottom quark pro¬ 
duction cross section in hadron reactions described above. The analysis of T(r 5) decays 
shows nevertheless that tighter experimental bounds on the bottom squark fraction are 
potentially powerful for the establishment of lower bounds on and nig. 

The hadronic width of the T is calculated in conventional QCD perturbation theory 
from the three-gluon decay subprocess, T —3g, and r 3 ^ a^. The SUSY subprocess 


adds a new term to the hadronic width from T 5 -f 5. If this new subprocess is present 
but ignored in the analysis of the hadronic width, the true value of = ^b) will be 

smaller than that extracted from a standard QCD fit by the factor (1 — FsusY/rsg) 3. 
For a contribution from the bb final state that is 25% of r^(lS) {i.e., a ratio rL = 10 
in Fig. 2), the value of extracted will be reduced by a factor of 0.9, at the lower 
edge of the approximately 10% uncertainty band on the commonly quoted value of 
as{mi,) p4l]. A thorough analysis would require the computation of next-to-leading 


order contributions in SUSY-QCD to both the 3g and bb amplitudes and the appropriate 
evolution of (Xs{iu) with inclusion of a light gluino and a light bottom squark. 

Direct observation of Upsilon decay into bottom squarks requires an understanding 
of the ways that bottom squarks may manifest themselves, discussed in more detail 
below. Possible baryon-number-violating R-parity-violating decays of the bottom squark 
lead to u + s; c + d; and c + s final states. These final states of four light quarks should 
be distinguishable from conventional hadronic final states mediated by the three-gluon 
intermediate state. If the b lives long enough, it will pick up a light quark and turn into a 
fi-mesino, B. Charged fi-mesino signatures in T decay include single back-to-back equal 
momentum tracks in the center-of-mass; measurably lower momentum than lepton pairs 
( < 4 GeV/c vs. ~ 5 GeV/c for muons and electrons); 1 -|-cos^0 angular distribution; 
and ionization, time-of-fiight, and Cherenkov signatures consistent with a particle whose 




FIGURE 2. Loci in the mj^-nig plane for which the rate for Upsilon decay into a pair of bottom squarks 
is either 10 times Lfor each T {nS), with n = 1 — 3 (lower set of curves) or equal to L(upper set). Above 
the curves, the rate would be less. Within each set of curves, the order is (bottom to top) 15,25,35. 


mass is heavier than that of a proton. At stake is discovery, or new limits on the mass, of 
the b as well as measurement of or new limits on the R-parity violation couplings of the 
b. 

Pseudoscalar decay into a pair of bottom squarks is forbidden, but the higher-order 
process of decay into a pair of fi-mesinos can proceed. Decays of the and into 
a pair of bottom squarks are allowed. 


Hadron Reactions 

Among the predictions of this SUSY scenario, the most clearcut is pair production of 
like-sign charged B mesons at hadron colliders, 5+5+ and B B . To verify the under¬ 
lying premise, that the cross section exceeds expectations of conventional perturbative 
QCD, a new measurement of the absolute rate for b production in run II of the Tevatron 
is important. A very precise measurement of % in run II is desirable. Since the fraction 
of b's from gluinos changes with pjb, a change of % is expected when the cut on pjb is 
changed. The b jet from g decay into bb will contain the b, implying unusual material 
associated with the b in some fraction of the bb data sample. The existence of light Vs 




means that they will be pair-produeed in partonic proeesses, leading to a slight inerease 
(~ 1 %) in the hadronie dijet rate. 

The SUSY approaeh inereases the b production rate at HERA and in yy collisions at 
LEP by a small amount, not enough perhaps if early experimental indications in these 
cases are confirmed [Q, Eull NEO SUSY-QCD studies should be undertaken. In these 
two cases, the apparent discrepancy may find at least part of its resolution in the fact that 
bb production occurs very near threshold where fixed-order QCD calculations are not 
obviously reliable. Uncertain parton densities of photons may play a significant role. 


Running of 

The presence of a light gluino and a light bottom squark slow the running of the strong 
coupling strength a^(;u). Above gluino threshold, the P function of (SUSY) QCD is 

a} / 2 1 \ 

P(o^5) = ^ f “11 + + 2 ] • (3) 

The b (color triplet scalar) contributes little to the running, equivalent to that of 1/4’th of 
a new flavor, but the g (color octet fermion) is much more significant, equivalent to that 
of 3 new flavors of quarks. A precise determination of ^( 0 X 5 ) appears to be the best way 
to exclude the presence of a light gluino or to establish its possible presence. 

In the standard model, a global fit to all observables provides an indirect measure¬ 
ment of QLs at the scale of the Z boson mass Mz- The value dsiMz) ~ 0.1184 ± 0.006 
describes most observables properly [^. Extrapolation from Mz to a lower scale /r with 
inclusion of a light gluino reduces <ls{b) from its pure QCD value. The presence of a 
light gluino, with or without a light bottom squark, also requires reanalysis of the phe¬ 
nomenological determinations of QLs{b) ^1 scales to take into account SUSY processes 
and SUSY-QCD corrections to the amplitudes that describe the relevant processes. To 
date, a systematic study of this type has not been undertaken, but, as mentioned above, 
consistency is achieved for T decays. A lesser value of as{mt,) leads, under slower evo¬ 
lution, to the same as{Mz). 

Slower running of (Xs{Q) also means a slower evolution of parton densities at small 
X, an effect that might be seen in HERA data for Q > mg. Presence of a scalar b in the 
proton breaks the Callan-Gross relation and yields a non-zero leading-twist longitudinal 
structure function Fl(x, Q) at leading-order. 


Z?-onia 

Bound states of bottom squark pairs could be seen as 7^ = 0+, 1 , 2+, ... mesonic 
resonances in yy reactions and in pp formation, with masses in the 4 to 10 GeV range. 
They could show up as narrow states in the invariant mass spectra at hadron 

colliders, between the //'P and T. At an collider, the intermediate photon requires 
production of a 7^^ =1 state. Bound states of low mass squarks with charge 2/3 were 


studied with a potential model pi|]. The small leptonic widths were found to preclude 
bounds for > 3 GeV. For bottom squarks with charge —1/3, the situation is more 
difficult. 


b lifetime and observability 


Strict R-parity conservation in the MSSM forbids b decay unless there is a lighter 
supersymmetric particle. R-parity-violating and lepton-number-violating decay of the b 
into at least one lepton is disfavored by the CLEO data [j^ and would imply the pres¬ 
ence of an extra lepton, albeit soft, in some fraction of Z? jets observed at hadron colliders. 
The baryon-number-violating R-parity-violating (/?p) term in the MSSM superpotential 
is U[ and are right-handed-quark singlet chiral superfields; and 

i,j,k are generation indices. The limits on individual and baryon-number violating 
couplings are relatively weak for third-generation squarks [^, < 0.5 to 1. 


The possible decay channels for the Z? are 123 : Z? 
223 '.b^c + s. The hadronic width is 


u + s', 213 '.b^c + d', and 


r(Z>-^jet+jet) = ^sin2 0~£|X"3|2. (4) 

j<k 


If = 3.5 GeV, T{b —> ij) = 0.08|?i"-3p GeV. Unless all are extremely small, the 

b will decay quickly and leave soft jets in the cone around the b. Z^-jets with an extra c 
are possibly disfavored by CDF, but a detailed simulation is needed. 

If the b is relatively stable, the b could pick up a light u or d and become a B or 5® 
“mesino" with J = 1/2, the superpartner of the B meson. The mass of the mesino would 
fall roughly in the range 3 to 7 GeV for the interval of b masses we consider. The charged 
mesino could fake a heavy muon if its hadronic cross section is small and if it survives 
passage through the hadron calorimeter and exits the muon chambers. Extra muon-like 
tracks would then appear in a fraction of the bb event sample, but tracks that leave some 
activity in the hadron calorimeter. The mesino has baryon number zero but acts like 
a heavy p - perhaps detectable with a time-of-fiight apparatus. A long-lived b is not 
excluded by conventional searches at hadron and lepton colliders, but an analysis [ ]27| ] 
similar to that for g’s should be done to verify that there are no additional constraints on 
the allowed range of b masses and lifetimes. 
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